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In  this  paper,  a  theoretical  model  for  the  performance  monitoring  and  fault  detection  of  fuel  ejectors  in 
the  hybrid  solid  oxide  fuel  cell  (SOFC)  system  is  proposed.  The  procedures  of  using  the  model  to  analyze 
ejector  properties  such  as  the  primary  mass  flow  rate,  the  secondary  mass  flow  rate,  the  recirculation 
ratio  and  steam  to  carbon  ratio  (STCR)  are  introduced.  Based  on  the  model,  the  anode  gas  recirculation 
performances  of  a  hybrid  SOFC  system  are  studied  under  various  operating  conditions.  Results  show  that 
the  model  can  be  used  to  evaluate  the  performance  of  ejector  not  only  in  the  critical  mode  but  also  in  the 
subcritical  and  back  flow  modes,  which  is  especially  useful  at  SOFC  off-design  operating  conditions  such 
as  start  up,  load  changes  and  shut  down. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells,  which  convert  chemical  energies,  such  as  hydrogen, 
methanol,  ethanol,  formic  acid,  and  methane,  into  electric  energy, 
possess  several  potential  advantages  over  conventional  combus¬ 
tion  power  generation  processes  including  higher  efficiency,  lower 
emissions,  and  higher  power  density.  Fuel  cell  technique  has  been 
considered  to  be  one  of  the  best  candidates  for  the  future  power 
generation  and  has  drawn  intensive  research  interests  in  recent 
years  [1,2].  Among  several  fuel  cell  schemes,  solid  oxide  fuel  cell 
(SOFC)  is  a  potential  alternative  in  the  distributed  power  genera¬ 
tion  for  domestic,  commercial  and  industrial  sectors.  In  the  anode 
side  of  an  SOFC  stack,  the  exhaust  is  rich  in  steam  (about  40-45%  in 
mass)  and  high  in  temperature  (around  900  °C),  which  can  be  recy¬ 
cled  by  means  of  an  ejector  to  provide  enough  steam  to  prevent 
carbon  deposition  and  sufficient  heat  for  endothermic  reforming 
reactions  in  the  cell  and  reformer  [3,4],  Another  ejector  can  also  be 
applied  in  the  cathode  side  to  recover  part  of  the  exhausted  gases 
to  replace  a  high  temperature  recuperator  [5,6].  In  these  ejectors, 
a  high  pressure  gas  called  as  the  primary  flow  is  used  to  entrain 
the  anode  or  cathode  exhausts.  Since  the  ejector  is  an  important 
component  in  the  recirculation  cycle  of  an  SOFC  system,  a  clear 
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understanding  of  its  working  principle  and  performance  character¬ 
istic  is  required. 

Ejectors  have  been  under  investigation  for  many  years,  and  sev¬ 
eral  modeling  techniques  are  available,  especially  in  the  areas  of 
refrigeration  and  chemical  engineering  [7-9],  Due  to  the  features 
of  high  recirculation  ratio,  low  pressure  increment  and  over  heated 
working  gases,  however,  very  few  existing  modeling  techniques  can 
be  directly  applied  to  SOFC  systems. 

Marsano  et  al.  [  10]  developed  an  SOFC  ejector  model  by  using  1 D 
modeling  technique.  The  developed  model  can  be  used  to  deal  with 
both  on-design  and  off-design  performance  evaluations.  However, 
this  model  does  not  account  for  the  operation  conditions  when  the 
amount  of  entrained  anode  recycle  gas  is  very  small  at  very  low 
primary  flow  pressure.  Later  on,  Ferrari  et  al.  [11]  improved  the 
modeling  technique  by  dividing  the  ejector  into  serials  of  calcula¬ 
tion  cells  where  the  governing  equations  were  numerically  solved 
by  1D-CFD  method  as  well  as  “lumped  volume”  treatment.  This 
improvement  makes  the  model  capable  of  predicting  the  transient 
behavior.  Recently,  Zhu  et  al.  [12]  proposed  a  new  modeling  tech¬ 
nique  for  fuel  ejectors  by  employing  a  2D  function  to  compute  fluid 
velocity  near  the  ejector  inner  walls.  The  developed  model  can 
be  applied  in  both  geometry  design  and  performance  simulation 
of  fuel  ejectors.  However,  the  model  cannot  accurate  predict  the 
ejector  performance  when  it  works  at  low  primary  flow  pressure. 

According  to  the  primary  flow  pressure,  the  ejector  perfor¬ 
mance  can  be  divided  into  three  operational  modes,  i.e.,  back  flow. 
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Nomenclature 

A  area(m2) 

Cp  specific  heat  of  gas  at  constant  pressure  (J  kg-1  K-1 ) 

D  diameter  (m) 

FC  fuel  cell 

k  specific  heat  ratio  of  gas 

m  mass  flow  rate  (kgs-1) 

M  Mach  number 

Mmj  Mach  number  of  the  mixing  layer  at  Section  3 

Mo  molecular  weight  (kg  mol-1 ) 

n  molar  flow  rate  (mol  s-1) 

nv  exponent  of  the  2D  curve 

P  pressure  (Pa) 

r,  R  radius  (m) 

Rs  gas  constant  (J  kg-1  K-1 ) 

Ru  universal  gas  constant  (J  mol-1 1<-1 ) 

T  temperature  (K) 

v,  V  velocity  (m  s-1 ) 

Greek  letters 

£exp  coefficient  accounting  for  friction  loss  during  the 
mixing  process 
p  density  (kg m-3) 

V p  isentropic  coefficient  of  primary  flow 

co  recirculation  ratio  {ms Imp) 

Subscripts 

P  primary  flow  (i.e.  inlet  fuel) 

S  secondary  flow  (i.e.  anode  recycling  gas) 

t  nozzle  throat 

0  ejector  inlet 

1  primary  flow  at  nozzle  throat 

2  nozzle  exit 

3  mixing  chamber  inlet 

4  mixing  chamber  outlet 

5  ejector  exit 

Superscripts 

i  chemical  component 


subcritical  and  critical  modes  [8,13].  The  ejector  may  work  in  the 
subcritical  mode  or  even  back  flow  mode  during  start  up,  load 
changes  and  shut  down.  In  these  cases,  the  ejector  performance 
characteristic  is  more  complex  than  that  in  the  critical  mode.  Noting 
that  both  the  heat  required  for  the  reforming  reactions  and  steam 
for  avoiding  carbon  deposition  are  supplied  by  the  entrained  anode 
gas,  unexpected  fluctuations  in  fuel  cell  system  could  occur  and 
then  the  system  might  run  under  some  “dangerous”  conditions. 
Therefore,  it  is  essential  to  develop  a  simple  and  accurate  fuel  ejec¬ 
tor  model  for  the  performance  monitoring  and  fault  detection  of 
fuel  ejectors  in  all  the  three  operational  modes. 

In  this  paper,  we  aim  to  develop  a  theoretical  model  for  fuel 
ejectors  in  the  hybrid  SOFC  system  for  performance  monitoring  and 
fault  detection  in  all  the  three  operational  modes.  Governing  equa¬ 
tions  for  computing  the  mass  flow  rate,  recirculation  ratio  and  STCR 
are  first  derived  based  on  the  thermodynamic  and  fluid  dynamic 
principles.  A  method  to  determine  the  two  key  parameters  PPE  and 
Ppc  is  proposed.  The  model  applications  in  the  performance  mon¬ 
itoring  and  fault  detection  for  fuel  ejectors  are  discussed.  A  hybrid 
SOFC  system  integrated  with  a  fuel  ejector  at  the  anode  side  is 
also  studied  using  the  proposed  model.  The  anode  gas  recircula¬ 
tion  behaviors  in  all  the  three  operational  modes  are  obtained  and 


analyzed.  Results  show  that  the  model  can  be  used  to  evaluate  the 
performance  of  ejectors  not  only  in  the  critical  mode  but  also  in  the 
subcritical  and  back  flow  modes.  This  property  is  especially  useful 
to  analyze  the  performance  of  an  off-design  operating  SOFC  system. 


2.  Fuel  ejector  in  anode  gas  recirculation  SOFC  system 


2.1.  SOFC  system  description 


A  typical  anode  gas  recirculation  SOFC  system,  which  mainly 
consists  of  three  components:  an  ejector,  a  reformer  and  a  fuel  cell 
stack,  is  schematically  shown  in  Fig.  1.  High  pressure  fuel  (primary 
flow)  passes  through  the  ejector  to  entrain  the  low  pressure  anode 
exhaust  (secondary  flow).  The  primary  flow  and  the  secondary  flow 
mix  in  the  mixing  chamber.  The  mixed  stream  shocks  in  the  diffuser 
and  then  enters  into  the  connected  reformer.  Inside  the  reformer, 
highly  endothermic  reactions  take  place: 

Reforming  :  CH4  +  H20  CO  +3H2  (fa) 

Shifting  :  CO  +  H20  C02  +  H2  (lb) 

The  reformed  fuel  is  fed  to  the  anode  side  of  FC  stack,  while  air 
is  supplied  to  the  cathode  side.  In  the  cathode,  oxygen  ions  passing 
through  the  electrolyte  layer  react  with  hydrogen,  and  the  electrons 
are  released.  These  electrons  pass  through  the  external  circuit  and 
reach  the  cathode  electrolyte  layer  to  make  the  circuit  close.  The 
reactions  inside  the  FC  stack  can  be  summarized  as  follows: 


Reforming  :  CH4  +  H20  CO  +3H2  (2a) 

Shifting  :  CO  +  H20  C02  +  H2  (2b) 

Electrochemical :  H2  +  |02  ->  H20  (2c) 


In  the  anode  gas  recirculation  SOFC  system,  STCR,  which  is  a  very 
important  parameter  to  evaluate  carbon  deposition  in  the  reformer 
and  FC  stack,  is  defined  as  follows  [10] 


(3) 


In  terms  of  fuel  ejectors  in  the  anode  gas  recirculation  cycle,  the 
STCR  can  be  computed  by: 


njff/yte  oMo' 

(nco/^uo^j  +  ^/^n'poMo^ 
where  the  recirculation  ratio,  co,  is  defined  as 


where  m$,  mp  are  the  mass  flow  rates  of  the  primary  flow  and  the 
secondary  flow,  respectively. 


2.2.  Ejector  operational  modes 


The  ejector’s  recirculation  ratio  is  strongly  influenced  by  three 
pressures:  primary  flow  pressure,  secondary  flow  pressure  and 
back  pressure  (pressure  of  gas  in  the  reformer).  Since  the  system 
load  of  an  SOFC  is  usually  adjusted  through  the  primary  flow  pres¬ 
sure,  its  effect  on  the  recirculation  ratio  is  shown  in  Fig.  2  [13,14]. 
Accordingly,  the  ejector  performance  can  be  divided  into  three 
operational  modes:  back  flow,  subcritical  and  critical.  The  primary 
mass  flow  rate  increases  with  the  primary  flow  pressure  in  all  the 
three  modes.  In  contrast,  the  behavior  of  secondary  flow  is  different 
in  each  mode: 
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Backflow  mode:  At  low  primary  flow  pressure,  no  secondary  flow 
is  entrained  into  the  ejector. 

Subcritical  mode:  The  ejector  starts  entraining  the  secondary  flow 
as  the  primary  flow  pressure  rises  to  Ppp.  The  recirculation  ratio  is 
very  sensitive  to  the  primary  flow  pressure  in  the  subcritical  mode. 
Critical  mode:  The  secondary  flow  shocks  in  the  ejector  and  reaches 
the  critical  mode  when  the  primary  flow  pressure  is  equal  to  Ppc. 
In  the  critical  mode,  the  secondary  flow  rate  decreases  first  and 
then  is  near  constant  in  the  high  pressure  region. 


her  is  lower  than  that  of  the  secondary  flow.  There  is  only  one  shock 
in  the  diffuser  (single-choking  [8]). 

In  the  critical  mode  (PP0  >  PPC)  as  described  in  Fig.  3(c),  the  pri¬ 
mary  flow  expands  after  the  nozzle  exit  introducing  a  series  of 
oblique  shocks  in  the  suction  chamber,  and  accelerates  the  sec¬ 
ondary  flow  to  choking  condition  at  the  mixing  chamber,  then  the 
mixed  flow  shocks  again  in  the  diffuser.  This  phenomenon  is  known 
as  double-choking  [8].  Since  the  secondary  flow  shocks  at  the  mix¬ 
ing  chamber  inlet,  the  ejector  working  at  the  critical  mode  is  more 
stable  than  the  other  two  modes. 


Detailed  flow  field  and  pressure  distribution  of  an  ejector  in  the 
three  operational  modes  are  shown  in  Fig.  3(a)-(c),  respectively. 
In  Fig.  3(a),  the  ejector  works  at  the  back  flow  mode  as  the  pri¬ 
mary  flow  pressure  ranges  from  0  to  PPg.  The  pressure  in  the  mixing 
chamber  is  higher  than  that  of  the  secondary  flow,  resulting  in 
part  of  the  primary  flow  is  reversed  and  no  shock  occurs  in  the 
ejector. 

In  the  subcritical  mode  (PPc  >PP,o  >Ppe).  the  secondary  flow  is 
entrained  into  the  ejector  due  to  the  pressure  in  the  mixing  cham- 


3.  Theoretical  model  development 

In  this  study,  it  is  assumed  that  the  ejector  meets  the  following 
conditions: 


1.  Both  the  primary  and  the  secondary  flows  are  ideal  gas  inside 
adiabatic  ejector  walls. 

2.  The  isentropic  relations  are  used  for  simplicity  in  deriving  the 
model. 

3.  The  primary  flow  velocity  is  uniform  in  the  radial  direction,  while 
velocity  of  the  secondary  flow  is  non-uniformly  distributed 
inside  the  ejector. 

4.  The  primary  flow  is  fully  heated  to  the  temperature  of  the 
secondary  flow  and  the  lost  heat  energy  of  secondary  flow  is 
negligible  (i.e.  rs,3  =  Ts,0;  TP,3  =  Ts,0)  [12], 

5.  Pressure  and  temperature  of  both  the  primary  and  the  secondary 
flows  are  uniformly  distributed  in  the  radial  direction  of  the  ejec¬ 
tor. 

6.  The  frictional  loss  during  the  mixing  process  is  taken  into 
account  by  a  coefficient. 


3.1.  Relations  between  ejector  inlet  and  Section  3 


Using  the  isentropic  flow  laws  and  the  energy  balance  equa¬ 
tion,  the  mass  flow  rate  of  the  primary  flow  through  the  nozzle 
is  [12] 


mP  -  AtpP;o(Vrpkp%pl'p,o)0'5 


V  kp  +  1 , 


(kp+l)/(2(kp-l)) 


(6a) 
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Rg,P  = 


Ru  XXo 
EnUMoi 


^nj,  0Mo' 

RuEnp.o 


Rg,pTp,o  Tp0 


for  a  mixture  inlet  fuel,  respectively. 

Relations  between  the  Mach  number,  velocity,  flow  diameter  of 
the  primary  flow  at  Section  3  and  the  inlet  boundary  conditions 
Pp,o,  Ps,o.  7s, o  can  be  written  as 


(Pp1 o 

U,o. 


MP’3  =  ^ 


(fcp-i)/fcp 


Op, 3  = 


V  fcp 

DtMp,3-°-5  f  2  +  (kP-  1)Mp,3: 


fexp 


2  +  (kP  -  1) 


2  \  (fep+l)/4^-!) 


(7) 


(9) 


where  Mp3,  and  Vp3  are  the  Mach  number  and  velocity  of  the  pri¬ 
mary  flow  at  Section  3,  respectively;  DP,3  is  the  diameter  of  the  flow 
area  of  the  primary  flow  at  Section  3.  The  details  of  deriving  Eqs. 
(7)  and  (9)  are  presented  in  Appendix  B  of  Ref.  [12], 

In  Section  3,  the  primary  flow  and  the  secondary  flow  are  sep¬ 
arated  by  a  mixing  layer  [12],  Inside  the  layer  is  the  primary  flow 
which  is  assumed  to  have  a  constant  velocity  in  the  radial  direction, 
and  the  secondary  flow  is  outside  the  layer  with  nonlinear  veloc¬ 
ity  distribution.  A  2D  velocity  function  for  the  primary  flow  and 
secondary  flow  in  Section  3  is  defined  as  follows  [12]: 


vr 


_  f  Vp  3  (0  <  r  <  RP  3) 

-  \w!3(l-r/R3)1/nv  (PP,3  <r<P3) 


(10) 


where  Rp,3  means  the  radius  of  the  mixing  layer  in  Section  3;  R3  is 
the  radius  of  the  mixing  chamber;  nv  is  the  exponent  of  the  veloc¬ 
ity  function.  Considering  the  velocity  and  the  radius  of  the  mixing 
layer  are  vr  =  Mml ^/kRgTS  0  and  r  =  RP3  =DP3/ 2,  respectively  [12], 
and  substituting  these  values  and  1/P,3  =  /WP  3  ^//<RgTs,o  into  Eq. 
(10),  we  obtain 


nv  = 


ln(l  -RP,3/R3) 
ln(Mml/Mp,3) 


(11) 


By  expressing  the  mean  mass  flow  rate  of  the  secondary  flow  at 
Section  3  as 


.-f 

■JRp  - 


pvrAA 


modes,  (a)  Back  flow,  (b)  subcritical  mode  and  (c)  critical  mode. 


where  i/tp  is  the  coefficient  relating  to  the  isentropic  efficiency  of 
the  primary  flow;  Tp,o  is  the  fuel  inlet  temperature;  the  average  gas 
constant  and  density  of  the  inlet  fuel  are  defined  as 


it  then  can  be  obtained  through  evaluating  the  integral  of  Eq.  (12), 

ms  _  (l  - 

where  the  average  density  of  the  secondary  flow  is  given  by 
nl  nMo! 

—  Ps,o  Ps,o  i 
Ps’°  ~  Rg,sTs,o  ~  7S, 0  RuY^nlo 


(13b) 
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3.2.  Mixing  process  between  Sections  3  and  4 

The  primary  flow  mixes  with  the  secondary  flow  in  the  mixing 
chamber.  The  mass,  energy  and  momentum  conservation  equations 
for  the  mixing  process  are  introduced  as  follows: 


Mass  conservation  equation: 


P4V4A3 

^g,m^4 


mP  +  ms 


(14) 


where  P4,  T4  and  V4  are  the  pressure,  temperature  and  velocity  of 
the  mixed  flow  at  Section  4,  respectively:  the  average  gas  constant 
of  the  mixed  flow  is 


*-£< 

^n^Mo* 


(15) 


Energy  conservation  equation: 

( mp  +  ms)(CpT4  +  |  V|)  =  mpCpTp  0  +  msQ>fs,o  (16) 

Momentum  conservation  equation: 

P4A3  +  ~(mP  +  ms)  =  {m  pVp,3  +  Pp,3^p,3)  +  (rns  Vs,3  +f*s,3^s,3) 

(17) 


where  the  specific  heat  ratio  of  the  fuel  in  the  mixing  chamber  is 
defined  as 


Note  that  some  geometry  properties  of  the  ejector  such  as  the 
converging  and  diverging  angles  and  the  chamber  lengths  are  not 
accounted  in  the  model  due  to  the  employment  of  ID  modeling 
technique.  Therefore,  the  influence  of  these  properties  on  the  ejec¬ 
tor  performance  is  not  revealed  in  this  model.  Fortunately,  this 
limitation  is  not  critical  since  these  geometry  properties  will  not 
seriously  affect  on  the  ejector  performance. 

4.  Operational  mode  estimation 

The  mass  flow  rate  of  the  primary  flow  and  secondary  flow  can 
be  analyzed  by  the  set  of  equations  from  Eqs.  (6)  to  (22).  Substi¬ 
tuting  these  values  into  Eqs.  (4)  and  (5),  a  theoretical  fuel  ejector 
model  for  the  anode  gas  recirculation  performance  evaluation  is 
finally  constructed.  As  stated  in  Fig.  2,  the  ejector  performance  is 
divided  into  the  three  operational  modes  by  the  two  values  of  Ppg 
and  PPC;  and  the  anode  gas  recirculation  behavior  is  quite  differ¬ 
ent  among  the  three  operational  modes.  In  order  to  better  monitor 
the  fuel  ejector  performance,  it  is  necessary  to  estimate  first  which 
mode  the  fuel  ejector  is  in,  and  then  to  evaluate  the  performance 
parameters  such  as  the  recirculation  ratio  and  STCR. 


where  <pm  expresses  a  coefficient  taking  account  the  frictional  loss 
in  the  mixing  process;  As,3  is  the  flow  area  of  the  secondary  flow  at 
Section  3  (As,3  =A3  -AP>3);  VS3  and  Ps, 3  which  are  the  velocity  and 
pressure  of  the  secondary  flow  in  Section  3,  can  be  determined  by 
Eqs.  (18)  and  (19),  respectively: 


Considering  Vs,3  =  ttis/(ps,o^s,3).  and  invoking  Eq.  (13a),  we 
have 


using  the  isentropic  flow  and  energy  conservation  law  for  the  sec¬ 
ondary  flow  from  Sections  0  to  3,  we  have 


3.3.  Pressure  diffusing  from  Section  4  to  Section  5 


In  the  subcritical  and  critical  modes,  the  mixed  flow  will  shock 
at  the  end  of  the  mixing  chamber  and  in  the  diffuser.  After  this 
shock,  the  kinetic  energy  of  the  mixed  flow  is  converted  into  pres¬ 
sure.  Assuming  the  pressure  diffusing  is  an  isentropic  process,  the 
pressure  at  the  exit  of  the  diffuser  P5  can  be  expressed  by 


Ps 

Pa 


km  —  1 
2 


(20) 


The  Mach  number  of  the  mixed  flow  at  Section  4  is  calculated  using 
the  following  relation: 


\Am%m T4 


(21) 


4.1.  Determine  Ppe 


The  primary  flow  starts  to  entrain  the  secondary  flow  when 
the  primary  flow  pressure  reaches  Ppp.  In  this  condition,  the  mass 
flow  rate  of  the  secondary  flow  is  zero,  i.e.  ms  =0;  Vs>3  =  0;  Ms, 3  =  0. 
Substitute  these  values  into  Eqs.  (14)  and  (16)-(19),  the  govern¬ 
ing  equations  for  the  flow  in  the  mixing  chamber  are  updated  as 
follows: 


P4V4A3  _ 
Pg,mT4 


(23) 


mp(Cpr4  +  \Vj)  =  mpCpTpo  (24) 

P4A3  -h'  =  (mpVpi3  +  Ps,oAP;3)  +  (Ps,t>As,3)  (25) 

Ps,3=Ps,o  (26) 


Vs, 3  =  0  (27) 

Due  to  the  non-linearity  and  close  coupling  of  these  equations, 
an  iterative  procedure  is  required  for  determining  Ppp  as  shown  in 
Fig.  4(a). 


4.2.  Determine  Ppc 


The  fuel  ejector  works  in  the  critical  mode  as  the  primary  flow 
pressure  is  greater  than  PPC.  In  this  operating  condition,  the  sec¬ 
ondary  flow  is  accelerated  by  the  primary  flow  and  always  shocks 
at  the  mixing  chamber  inlet.  At  the  critical  mode,  it  can  be  reason¬ 
ably  assumed  that  only  the  layer  between  the  primary  flow  and 
secondary  flow  in  Section  3  is  in  the  choking  condition,  i.e.  Mm]  =  1 
[12],  Substituting  this  into  Eq.  (11),  we  have 


ln(l-Rp>3/R3) 

ln(l/MP,3) 


(28) 


Using  this  condition,  a  detailed  calculation  flowchart  for  deter¬ 
mining  Ppc  is  given  in  Fig.  4(b). 
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5.  Application  procedures 

For  a  given  ejector  geometry,  the  anode  gas  recirculation  per¬ 
formance  depends  on  Ppo  and  Ps,o.  Pp.o  and  Ts  o,  and  the  chemical 
composition  of  inlet  fuel  and  anode  recycle  gas.  The  model  appli¬ 
cation  procedure  in  performance  monitoring  and  fault  detection  is 
introduced  as  follows: 

(1) .  Estimate  which  operational  mode  the  ejector  is  in:  the  ejector 

performance  is  divided  into  the  three  operational  modes  by 
PPE  and  PpC,  which  can  be  determined  by  the  flowchart  shown 
in  Fig.  4.  Operational  mode  of  the  fuel  ejector  is 

Operational  mode 

back  flow  mode;  (0  <  Pp,o  <  Ppe) 
subcritical  mode;  (PPE  <  PP  0  <  Ppe)  (29) 

critical  mode;  (Pp,o  >  Ppe) 

(2) .  Compute  mP  and  m $:  the  mass  flow  rate  of  the  primary  flow 

mP  can  be  computed  from  Eq.  (6a)  as  the  fuel  ejector  works 
in  all  the  three  modes.  But  for  computing  the  mass  flow  rate 
of  the  secondary  flow  ms,  the  governing  equation  is  different 
in  the  three  modes.  In  addition,  as  shown  in  Fig.  2,  it  is  found 
that  ms  is  approximately  linear  with  the  primary  flow  pressure 
in  the  subcritical  mode  and  satisfies  the  boundary  conditions; 
Ppo  =  PpE,  m$  =  0  and  Ppo  =  Ppc,  m$  =/(Ppc).  In  order  to  reduce 
the  calculation  time,  a  linear  function  is  defined  to  approach 
the  secondary  mass  flow  rate  in  the  subcritical  mode  by  con¬ 
sidering  the  boundary  conditions.  Governing  equations  for  the 
secondary  mass  flow  rate  in  the  three  operational  modes  are 
expressed  by: 

{-  (0  <  PP,o  <  Ppe) 

/(Ppc)%c°:;;eE  (Ppe  <  PP,o  <  Ppe)  (30) 

/(Pp,o)  (Pp,0  >  Ppe) 


where  /  is  the  function  of  m$.  Inputs  of  the  function  are  the 
ejector  geometries  Dt  and  D3 ,  and  the  operation  conditions  Ppo, 
Ps,o.  Tp,o  and  Tso;  and  output  is  ms.  The  detailed  computation 
procedure  is  given  in  Fig.  5. 

(3).  Determine  the  recirculation  ratio  and  STCR:  the  recirculation 
ratio  and  STCR  can  be  computed  from  Eqs.  (4)  to  (5)  once  mP 
and  ms  are  determined. 

6.  Results  and  discussions 

As  an  example,  Table  1  shows  the  typical  values  for  an  anode 
gas  recirculation  SOFC  system  analyzed  in  [12].  Starting  with  these 


Fig.  5.  Calculation  flowchart  for  ejector  performance  simulation. 


1128 


Y.  Zhu  etal.  /Journal  of  Power  Sources  185  (2008)  1122-1130 


Table  1 

On-design  values  of  the  ejector  and  SOFC 


Parameter  Value 


Fuel  inlet 

Composition  (molar,  %) 

CH4  100 

Flow  rate  (kgs-1)  0.0094 

Pressure  (bar)  10.06 

Temperature  (K)  673 

Anode  gas  recirculation  gas 
Composition  (molar,  %) 

H2  4.895 

CO  3.785 

H20  61.74 

C02  29.58 

Pressure  (bar)  3.8 

Temperature  (K)  1280 

Fuel  cell  parameters 

Cathode  thickness  (cm)  0.035 

Electrolyte  thickness  (cm)  0.017 

Anode  thickness  (cm)  0.030 

Overall  cell  area  (m2)  95 

Fuel  cell  operation  conditions 

Fuel  utilization  0.85 

FC  pressure  (bar)  3.80 

FC  pressure  loss  (kPa)  5.7 

Air  inlet  pressure  (bar)  3.84 

Air  inlet  temperature  (K)  1000 

Air  flow  rate  (kgs-1)  0.47 


data,  the  influences  of  the  fuel  inlet  conditions,  the  cell  operation 
pressure  and  temperature,  the  fuel  utilization,  etc.,  on  the  hybrid 
SOFC  system  are  carefully  studied  based  on  the  proposed  fuel  ejec¬ 
tor  model.  During  the  simulation,  the  SOFC  model  developed  by 
Costamagna  et  al.  [15]  is  adopted,  which  allows  the  evaluation  of 
both  on-design  and  off-design  behavior  of  the  SOFC  system.  The 
main  assumptions  in  this  study  are  listed  as  follow: 


1.  Temperature  within  all  the  components  of  SOFC  system  is  uni¬ 
formly  distributed. 

2.  Cathode  flow  is  composed  of  21%  O2  and  79%  N2.  Fuel  is  pure 
CH4. 

3.  The  reforming  and  shifting  reactions  are  at  equilibrium  in  the 
reformer  and  FC  stack. 

4.  Temperature  of  the  gases  at  the  outlet  of  the  reformer  and  FC 
stack  are  equal  to  the  reformer  and  FC  stack  temperature,  respec¬ 
tively. 

5.  The  pressure  loss  in  FC  stack  is  equal  to  1.5%  of  the  FC  operation 
pressure. 

6.  The  concentration  loss  is  fixed  at  1.5  x  10-7  £2. 


A  design  value  of  STCR  equals  to  2.4  and  must  be  higher  than  2  to 
avoid  the  carbon  deposition  [2,16],  The  fuel  utilization  coefficient 
Uf  is  define  as 


consumed 

Hz 


(31) 


where  n“nsumed  represents  the  reaction  rate  of  H2  in  the  FC  stack, 
and  n}^2 ,  ng,,  are  molar  flow  rate  of  H2,  CO  and  CH4  into  the 
FC  stack,  respectively. 


6.1.  Ppe  and  Ppe 


Ppe  and  PpC  are  affected  by  the  operating  conditions.  Effects  of 
the  primary  flow  temperature,  and  the  secondary  flow  tempera¬ 


ture  and  pressure  on  Ppe  and  PPC  are  presented  in  Fig.  6(a)-(c), 
respectively. 

Fig.  6(a)  shows  the  behavior  of  Ppe  and  PPC  at  different  primary 
flow  temperatures.  It  is  observed  that  both  Ppe  and  Ppe  increase  with 
the  primary  flow  temperature.  This  can  be  explained  as  follows: 


(1)  Higher  temperature  means  lower  density  for  a  fixed  pressure 
gas; 

(2)  Gas  dynamic  viscosity  increases  when  temperature  increases 
and  pressure  is  constant. 


Secondary  flow  temperature  (K) 


flow  temperature  and  (c)  secondary  flow  pressure. 
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Both  lead  to  less  mass  flow  rate  of  the  primary  flow  according 
to  the  fluid  dynamics  principles. 

The  relations  between  the  two  pressure  values  and  the  sec¬ 
ondary  flow  temperature  are  shown  in  Fig.  6(b).  Both  Ppe  and  Ppc 
decrease  with  the  secondary  flow  temperature,  which  implies  that 


Primary  flow  inlet  pressure  (bar) 


Primary  flow  inlet  pressure  (bar) 

Fig.  8.  Performance  of  a  hybrid  SOFC  at  on-design  conditions:  (a)  anode  gas  recir¬ 
culation  performance  and  (b)  fuel  cell  stack  performance. 


the  ejector  will  start  to  entrain  the  secondary  flow  earlier  if  the 
secondary  flow  temperature  is  higher. 

In  Fig.  6(c),  it  is  seen  that  Ppe  and  Ppc  increase  as  the  secondary 
flow  pressure  increases.  Noting  that  value  of  the  ejector  exit  pres¬ 
sure  (P5)  remains  1.015  times  of  the  secondary  flow  pressure,  we 
conclude  from  Fig.  6(a)  to  (c)  that  the  fuel  ejector  starts  to  entrain 
the  secondary  flow  once  the  primary  flow  pressure  is  1.11-1.19 
times  of  the  ejector  exit  pressure,  and  will  work  in  the  critical  mode 
as  the  primary  flow  pressure  is  greater  than  2.17-2.25  times  of  the 
ejector  exit  pressure. 

6.2.  Fuel  ejector  performance 

The  most  important  performance  parameters  in  the  fuel  ejector 
obviously  are  the  two  mass  flow  rates,  recirculation  ratio  and  STCR. 
With  the  new  model,  these  parameters  can  be  analyzed  for  all  the 
three  operational  modes. 

Keeping  the  primary  flow  temperature  constant  and  the  anode 
exhaust  condition  at  the  design  values  as  stated  in  Table  1,  the 
detailed  relationships  between  the  two  mass  flow  rates,  recircu¬ 
lation  ratio,  STCR  and  the  primary  flow  inlet  pressure  are  shown 
in  Fig.  7.  It  is  seen  that  the  mass  flow  rate  of  the  primary  flow 
increases  linearly  with  the  primary  flow  inlet  pressure  in  accor¬ 
dance  with  Eq.  (6a).  However  the  primary  flow  inlet  pressure  has 
distinct  influences  on  the  other  performance  parameters: 

•  When  the  primary  flow  pressure  (Ppe)  increases  up  to  about 
4.45  bar,  the  secondary  flow  starts  to  be  entrained  into  the  ejec¬ 
tor.  The  mass  flow  rate  of  the  secondary  flow  increases  gradually 
with  the  primary  flow  pressure  in  the  subcritical  mode,  so  do  the 
recirculation  ratio  and  STCR. 

•  The  ejector  works  in  the  critical  mode  when  the  primary  flow 
pressure  is  greater  than  8.47  bar.  The  mass  flow  rate  of  the  sec¬ 
ondary  flow  decreases  first  and  then  remains  quite  constant  in 
the  high  primary  flow  pressure  region. 

6.3.  Fuel  cell  performance 

First  keeping  the  fuel  cell  pressure  PFc  at  3.8  bar  and  the  fuel 
utilization  coefficient  Uf  at  0.85,  the  SOFC  system  performance 
is  studied  by  varying  the  other  operation  conditions  as  shown  in 
Fig.  8(a)  and  (b).  Then  varying  PFC  and  U{,  the  off-design  behaviors 
of  recirculation  ratio  and  STCR  are  shown  in  Fig.  9(a)  and  (b). 
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cell  operation  pressure  and  (b)  fuel  utilization  coefficient. 

Fig.  8(a)  shows  the  two  mass  flow  rates,  recirculation  ratio 
and  STCR  with  the  variation  of  the  primary  flow  inlet  pressure. 
These  performance  parameters  have  different  behaviors  but  similar 
trends  with  those  obtained  from  Fig.  7.  The  differences  are  caused 
by  the  temperature  and  chemical  composition  of  the  anode  recycle 
gas  which  are  determined  by  SOFC  operating  conditions  as  shown 
in  Fig.  8(a),  while  they  are  always  fixed  at  the  design  points  in  Fig.  7. 

In  Fig.  8(b),  the  performance  parameters  of  the  fuel  cell  stack: 
voltage,  power  and  current  density  are  illustrated.  It  is  found  that 
they  have  the  same  trends  in  both  the  subcritical  mode  and  the 
critical  mode  because  that  these  parameters  mainly  depend  on  the 
mass  flow  rate  of  the  inlet  fuel  while  Uf  keeps  at  the  design  point. 
The  current  density  is  proportional  to  the  fuel  inlet  pressure  in  the 
subcritical  and  critical  modes.  It  may  be  explained  that  the  flow  rate 
of  CH4  into  fuel  cell  stack,  njijj  ,  increases  with  the  fuel  inlet  pressure 
so  that  n“nsumed  will  increase  according  to  Eq.  (31 ),  resulting  in  an 
increase  of  the  current  density.  From  Fig.  8(b),  it  is  also  seen  that 
the  FC  power  increases  with  the  fuel  inlet  pressure  as  the  is 
increased  while  Uf  is  kept  constant.  In  all  the  simulated  fuel  inlet 
pressure,  FC  voltage  is  near  0.53  V. 

The  recirculation  ratio  and  STCR  with  the  variation  of  the  FC 
operation  pressure  is  shown  in  Fig.  9(a).  The  recirculation  ratio  and 


STCR  behaviors  in  the  critical  mode  differ  from  that  in  the  subcriti¬ 
cal  mode.  The  STCR  can  drop  below  the  limited  value  for  a  lower  FC 
pressure  in  the  critical  mode,  but  different  in  the  subcritical  mode 
due  to  the  FC  operation  pressure  has  an  effect  on  Ppe  and  Ppc  as 
shown  in  Fig.  6. 

Fig.  9(b)  shows  the  results  obtained  by  varying  the  fuel  utiliza¬ 
tion  coefficient  at  different  fuel  inlet  pressures.  The  figure  indicates 
that  the  STCR  is  strongly  influenced  by  Uf,  as  Uf  has  a  direct  effect 
on  the  conversion  rate  of  FI2,  which  affects  the  FC  performances 
such  as  temperature,  pressure  and  chemical  composition  at  the 
anode  recycle  gas.  Results  confirm  that  for  a  low  Uf,  the  STCR  can 
easily  drop  below  the  limited  value  and  thus  suffers  from  carbon 
deposition. 

7.  Conclusions 

A  theoretical  ejector  model  for  performance  monitoring  and 
fault  detection  in  the  hybrid  SOFC  system  was  proposed  in  this 
paper.  The  model  was  used  to  analyze  the  fuel  ejector  proper¬ 
ties  such  as  the  primary  mass  flow  rate,  the  secondary  mass  flow 
rate,  the  recirculation  ratio  and  STCR  not  only  in  the  critical  mode 
but  also  subcritical  and  back  flow  operational  modes.  Furthermore, 
a  method  for  determining  PPE  and  PPC,  consequently,  the  ejector 
operational  mode  was  introduced. 

By  utilizing  the  new  model,  the  performances  of  an  anode  gas 
recirculation  SOFC  system  integrated  with  a  fuel  ejector  were  inves¬ 
tigated  for  the  back  flow,  subcritical  and  critical  operational  modes. 
The  main  conclusions  from  the  simulation  are: 

1.  Two  parameters  Ppe  and  Ppc  which  divide  the  fuel  ejector  perfor¬ 
mance  into  three  operational  modes  should  be  updated  as  soon 
as  the  fuel  cell  operating  conditions  changes. 

2.  The  ejector  will  not  work  at  the  critical  mode  if  the  primary  flow 
inlet  pressure  is  less  than  PPC.  In  such  case,  recirculation  ratio 
and  STCR  need  to  be  carefully  monitored. 

3.  The  fuel  ejector  starts  to  entrain  the  secondary  flow  as  the 
primary  flow  pressure  is  about  1.15  times  of  the  ejector  exit  pres¬ 
sure:  reaches  the  critical  mode  when  the  primary  flow  pressure 
is  approximate  2.2  times  of  the  ejector  exit  pressure. 
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